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APPLICATION OF BOUNDARY INTEGRAL METHOD TO ELASTIC 
ANALYSIS OF V-NOTCHED BEAMS 
by Walter Rzasnicki, Alexander Mendelson, and Lynn U. Albers 

Lewis Research Center 

SUMMARY 

A semidirect boundary integral method, using Airy's stress function and its deriva- 
tives in Green’s boundary integral formula, is used to obtain an accurate numerical 
solution for elastic stress and strain fields in V-notched beams in pure bending. The 
proper choice of nodal spacing on the boundary is shown to be necessary to achieve an 
accurate stress field in the vicinity of the tip of the notch. Excellent agreement is ob- 
tained with the results of the collocation method of solution. 


INTRODUCTION 

Knowledge of the stress distribution in the neighborhood of a singularity, such as 
the tip of a V-notch in a beam loaded in tension or bending, is of fundamental importance 
in evaluating the resistance to fracture of structural materials. Elastic solutions to 
various geometries have been obtained by a number of different methods. Among the 
more effective ones, are the complex variable method (ref. 1), collocation method 
(ref. 2), and finite element method (ref. 3). However, the first two of these methods 
are not general enough nor readily adaptable to three-dimensional or elastoplastic prob- 
lems. And the finite element method requires solutions of large sets of equations and 
fails to give sufficiently fine resolution in the vicinity of crack tips. 

The recently developed boundary integral methods (ref. 4) offer an attractive alter- 
native to other methods of analysis. These methods have a number of advantages which 
may be listed as follows: 

(1) They obviate the need for conformal mapping. 

(2) Mixed boundary value problems are handled with ease. 

(3) Stresses and displacements are obtained directly without need for numerical dif- 
ferentiation. 



(4) No special considerations are needed for multiply connected regions. 

(5) The internal stresses and/or displacements are calculated only where and when 
needed. 

(6) The extension to three-dimensional problems is direct. 

(7) Nodal points are needed only on the boundary instead of throughout the interior 
as required by finite element methods. 

The last point, which is probably the most important one, is illustrated in figure 1. 
For the finite element method, the whole region must be covered by a grid producing a 
large number of nodal points and corresponding unknowns. Thus a large number of 
simultaneous equations must be solved. For the boundary integral methods, nodal 
points are taken only on the boundary, resulting in a much smaller number of unknowns. 

The object of this report is to present a solution to the elastic problem of a V- 
notched beam in pure bending using one of the boundary integral methods described in 
reference 4. A necessary part of this solution is the strategy by which the nodal spac- 
ing on the boundary is chosen to achieve stable and accurate solutions. It is intended to 
use the spacing which gives good elastic results in the extension of the method to the 
more interesting and complicated elastoplastic problem, which will be discussed in a 
subsequent report. 


SYMBOLS 


a 


c.., 1 
iy i y i] , ‘ 


d ij ,e ij ,f ij 

A ij, B ij> c ij> D ij? 

• E ij’ F ij’ G ij ,H ij’‘ 

V K ij 

C 
E 
J 

K x 


notch depth 

dimensionless notch depth, a/w 
coefficients in boundary equations 


coefficients in stress equations 

boundary contour 

Young’s modulus of elasticity 

value of Rice's contour integral 

stress intensity factor at notch tip for mode I crack opening 
half length of beam 
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n 

n, s 

P(x,y) 

q(C,T?) 

/V 

q 

R 


r ij 

r,B 

s 


T i 

u. 


u 


u(x, y), v(x, y) 
W(c) 


w 


x, y,z 
x,y 


a 


5 « 


e ij 

r 


M 

p 


a ij 

a 

max 

°x’ a y’ a z’ a xy 

a 0 

$ 

<P 


order of stress singularity in vicinity of tip of notch 

unit vectors normal and tangent to contour C 

point on contour C or in region R 

point on contour C 

dimensionless load, a max / a Q 

planar region bounded by closed contour C 

distance between two points having coordinates (x, y)^ and (^ , tj)j 

polar coordinate directions 

length measured along contour C 

stress vector active along boundary 

displacement vector 

displacement in y-direction 

arbitrary continuous functions 

strain energy density 

width of beam 

rectangular Cartesian coordinate directions 

dimensionless rectangular Cartesian coordinate directions, x/w, y/w 
notch angle, deg 
Kronecker delta 
strain tensor 

arbitrary contour surrounding crack tip 
Poisson's ratio 

rectangular Cartesian coordinate directions 
o 

function of r, p = r In r 
stress tensor 

maximum nominal bending stress 

components of stress tensor in Cartesian coordinates 

tensile yield stress 

2 

function of Airy stress function, $ = V cp 
Airy stress function 


V 2 Laplace's operator, 3 2 /3x 2 + 3 2 /3y 2 

V 4 biharmonic operator, 3 4 /3x 4 + 2 3 4 /3x 2 3y 2 + 3 4 /3y 4 

Subscripts : 

i,j,k integers 

Superscripts: 

~ dimensionless quantity 

' derivative in outward normal direction, 3/3n 


ANALYSIS 


, . The Integral Equation Method 

The problem of determining the state of stress in a plane elastostatic problem can 
be reduced to solving a homogeneous biharmonic equation 


V 4 <p(x,y) = 0 (1) 

subject to appropriate boundary conditions, where cp(x, y) is the Airy stress function 
and the components of the stress tensor, in Cartesian coordinates, are 


°x = 


_ 9 2 <P 


3y 


d 2 (p 


3 2 <p 


xy 


3x3y 


( 2 ) 


For the problem under consideration (fig. 2) the stress function <p (x, y) and its 
outward normal derivative 3 <^/3n must satisfy the following boundary conditions 
(ref. 5): 


4 


along boundary OA and OA’: 


'x 


along boundary AB and A'B': 


<P (x, y) = 0; l£ = 0 
Sn 


<p(x, y) = 0;*£ = 0 

3n 


along boundary BC and B'C': 


y ( 3 ) 


*<*, y) - ♦ a 2 x + 2!) + „ + « + A l£ = o 

w \ 3 3/ max y 2 2/ 3n 


along boundary CD and C'D: 


cr w 


V (x,y)=-iSH_;l5» = o 

6 3n 




The method of solution utilizes Green's second theorem to reduce the homogeneous 
biharmonic equation (1) to two coupled, Fredholm type, integral equations, which must 
satisfy the specified boundary conditions. 

Green's second theorem states 


If 


(uV 2 v - vV 2 u) dx dy = 



(4) 


where the arbitrary functions u(x, y) and v(x, y) and their derivatives of first and sec- 
ond order are continuous in the simply connected region R, bounded by a sectionally 
smooth curve C. The notation employed is depicted in figure 3. 

If we let 



( 5 ) 


then it follows, that 



vV 


<p ) dx dy = 



— (V 2 v) - V 2 v + V 2 <p iX - v — iy 2 (p) ds 
3n 3n 3n 3n 


( 6 ) 
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Let us introduce the function 4>(x,y) such that 

$ = V 2 (p (7) 

Substituting equation (7) into equation (1) yields 

V 2 $(x,y) = 0 (8) 

Let r(x, y; ^,77) be the distance between any two points P(x, y) and q in the 

region R, as shown in figure 3, such that PcR+C and q c C. 

Substituting relations (5), (7), and (8) into Green's second theorem, equation (4), 

o 

choosing first, v = In r (the fundamental solution of V v = 0) and taking into account the 
singularity at r = 0 (ref. 5) results in 


2tt$(x 


,y>- /V 

Jr, l dn 


(In r) - — In r 
8n 


ds for P c R 


(9) 


and 


7r$(x, y) = f ["$ — (In r) - — In 
Jr, L 9n 9n 


ds for P c C 


( 10 ) 


Combining equation (6) with equations (1) and (7) and choosing this time v = p - 

2 

r In r, and taking into account the singularity at r = 0, we obtain (ref. 6) 


87r<p(x,y) = f \<p i- (V 2 p) - V 2 p + $ i£ - It p 
Jr L 3n 9n 0n 0n . 


ds for Pc R (11) 


and 


4n-<p(x 


, y) = f \<p l(v 2 p)-liv 2 p + ^-l£p 

jJ, L 3n 0 n 0 n 0 n . 


ds for P c C (12) 


Equation (11) would give us directly a solution to the biharmonic equation (1) pro- 
vided the functions <p(x, y), d(p(x, y)/0n, V 2 <p(x, y), and 3 [v 2 ^(x, y)]/0n were known 
on the, boundary C. 

However, only the stress function <p and its outward normal derivative d<p/dn 


6 



2 2 

are specified on the boundary. The values of V cp = $ and 3(V cp)/dn = d$/dn on 
the boundary must be compatible with these boundary conditions. To assure this com- 
patibility, we have to solve the system of coupled integral equations (10) and (12), which 
contain the unknown functions <f> and d$/dn. 

Once the values of <f> and 9$/ 3n on the boundary C of the region R are known we 
can proceed with the calculation of the stress field in the region R utilizing equations (11) 
and (2). 


Numerical Procedures 

Solution of the integral equations . - Since it is generally impossible to solve the 
system of coupled integral equations analytically, a numerical method is utilized in 
which the integral equations (10) and (12) are replaced by a system of simultaneous al- 
gebraic equations. 

For simplicity of notation, let us denote the normal derivatives by prime super- 
scripts. 

Let us assume that the function <i> and its normal derivative <£' are piece- wise 
constant on the boundary C. We can divide the boundary into n intervals, not neces- 
sarily equal, numbered consecutively in the direction of increasing s. The center of 
each interval is designated as a node and assigned the constant values of and <3>J. 
The arrangements of boundary subdivisions is shown in figure 4. 

Using these assumptions, equations (10) and (12) can be replaced by a system of 
2n simultaneous algebraic equations with 2n unknowns, that is, <i>j and 


n 


= > (a.4>. + b..$!) 

i L ' i] 3 ij r 
j=l 


u 

Atrip. = > (c..$. + d..$! + e..<p. + f..<pl) 

Y i L K u i x] j irj irr 

i =1 


(13) 


4 

where i = 1, 2, 3 . . . n. The coefficients appearing in the boundary equations (13) are 
defined by the following relations: 


7 



‘ij = / (ln r ij ) ' 


ds 


% =-/ tar ij ds 

c y - f p;, ds 
•'j 

d.. = - / P-- ds 
ij i] 

•'j 

e ii - [ (A))’ ds 

v- / v V s 

J i ^ 


(14) 


XL i.1. 

where integration is taken over the j interval, and r.. is the distance from i node 
th " 

to any point in the j interval. The normal derivatives in equation (14) are taken on 
the j th interval. 

The coefficients given by equations (14) can be evaluated by Simpson's rule for 
i * j, and analytically for i = j. For the boundary intervals which can be represented 
by straight lines an analytical solution is possible for all the coefficients. 

Equations (13) expressed in matrix form become 


h - V] 

hi' 


V' 


[0] 

[»]' 


TrJ 

n x n 

n x n 


n x 1 


n x n 

n x n 


n x 1 

w 

h] 

◄ 

[*j] 

> = - 

[ e ij - V'] 

U 


M * 

n x n 

n x n 


n x 1 

L J 


n x n 

n x n 


n x 1 

L J 


(15) 
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Thus, the problem is reduced to the solution of the following matrix system: 

[B] {X} = {R} 


(16) 



where [B] is a 2n x 2n matrix and { x } and {r} are 2n x 1 column matrices. 

To calculate stresses for any point in the region R from the stress function (11), 
we need not perform any numerical differentiation of the stress function. Once <3? and 
4>’ are known on the boundary we can differentiate under the integral sign in this equa- 
tion and then obtain stresses by the same type of numerical integration as in equa- 
tions (13). Applying equations (2) to equation (11) and using notations of equations (14) 
yields the stress equations 


n 


''v 


8 ^ x (x,y)i = £ (A^j + + Cyttj + D^j) 

J=1 


8 ^ y (x,y)i = £ - B.^j + By*, + F..$j) 

1=1 


(17) 


n 


-8^ xy (x,y)i = £ <V) + H ij^j + *tj*J + VP 

j=l 


where i refers to any point in the stress field, and the coefficients A.j, B.^, C.j, D.j, 
Ejj, F^j, G.j, fbj, I.j, and are obtained by appropriate differentiation under the 
integral sign of the coefficients given by equations (14), and are listed in the appendix. 

The stress function cp is not constant on the loaded boundaries BC and B’C’. 

The assumption that it is piece-wise constant may lead to appreciable errors in the 
numerical results. To overcome this difficulty, the summations given in equations (13) 
and (17) for intervals lying on the loaded boundaries and involving the stress function are 
replaced by direct integration. 

Boundary interval size . - The number of nodal points prescribed for the boundary 
is theoretically unlimited. However, computer storage capacity for the computer used 
and the difficulties associated with inversion of large matrices limited the order of the 
coefficient matrix [B] of equation (16) used herein to 140. 

Because of geometric and loading symmetry about the x axis it is possible to re- 
duce the number of unknowns. For 2n total number of nodal points the number of 
equations and unknowns, and <f>?, is reduced from 4n to 2n. Additional reduction 
in the number of unknowns is accomplished by taking into consideration the St. Venant’s 
effect at the loaded boundaries. By definition 


4>. = V 2 <p. = a 
1 x- 


+ a 
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Since, near the loaded boundaries BC and B'C', a can be assumed to be essen- 

j 

tially zero if L ^ 1.2 (ref. 2), then it follows that 



<f> 


j 



an 


(18) 


Thus on the boundaries BC and B'C', and <f>j are known. 

The arrangement of boundary subdivisions and nodal points is shown in figure 5. 
Note that the corner points are always designated as interval points, never as nodal 
points, thus eliminating discontinuous functions from the numerical analysis. 

Since the vicinity of the crack tip is of greatest interest, a fine nodal spacing along 
the notch was chosen. To reduce the error introduced by the change in the interval size 
(ref. 7) around boundary points A and A' and at the same time to obtain fine resolu- 
tion at the tip of the notch, the boundary along the notch was divided into a number of 
intervals progressively decreasing in length. The rate of change in the interval length 
and the resulting length of the smallest interval was found to have a great influence on 
the stress field in the vicinity of the tip of the notch. 

Each of the boundaries AB, A'B', DC and DC' was divided into 5 intervals of 
equal length. Each of the boundaries BC and B'C' were assigned values of <f>j and 
<f>j, given by equations (18), at 15 uniformly spaced nodal points. To complete the 
boundary nodal arrangement, 60 nodal points were taken along each of the boundaries 
along the notch. The rate of change in the interval's length along these boundaries was 
optimized for all cases by the method illustrated in figure 6 for the case of a specimen 
with a 10° edge notch and a/w = 0.3. The dimensionless y-directional stresses at 
three locations below the tip of the notch were plotted as a function of tapering ratio, 
that is, the ratio of the lengths of two consecutive boundary intervals. Note that the 
stresses far enough from the tip were insensitive to changes of the taper. The general 
pattern was that the stress converges smoothly to its correct value as the tapering ratio 
increases until at a certain point it starts to oscillate erratically. This behavior is 
caused by the violent oscillations in $ and <f>' near the tip of the notch that occur when 
the minimum interval length becomes too small. 

This result suggests the following strategy in choosing the optimum tapering ratio 
for a given number of boundary intervals along the notch: 

(1) Plot the x-directional or y-directional stresses at several locations below the 
tip of the notch, where good resolution of stresses is desired, for a sequence of taper- 
ing ratios. 
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(2) Start with a ratio of 1. 00 and increase it uniformly. Continue plotting these 
curves from the smooth monotonic regime into the region where the variation becomes 
erratic. Pick as an optimum ratio the largest one that is in the smooth regime for all 
curves. 

For the case of a 10° edge notch and an a/w = 0. 3, an optimum tapering ratio of 
1. 09 is associated with the smallest dimensionless interval length of approximately 
0. 00015. For the other cases considered, optimum ratios were found to be in the range 
of 1. 08 to 1. 10. The resulting smallest dimensionless boundary interval length varied 
from 0. 0001 to 0. 0002. 

The nodal arrangement shown in figure 5 was used for all cases considered, result- 
ing in a set of 140 equations containing 140 unknowns, that is, <£. and i = 1, 2, 

3, . . ., 70. 

The numerical calculations were performed on an IBM 7094 digital computer using 
single -precision arithmetic. The matrix system given by equation (16) was solved using 
the modified Gauss elimination method, which utilizes pivoting and forward and back- 
ward substitutions. 


RESULTS AND DISCUSSION 

The computations were performed for cases of 10°, 30°, and 60° notch angles and 
varying notch depths. Selected results were also obtained for a 3° notch angle. Fig- 
ures 7 to 9 show x-directional and y-directional stress distribution as a function of dis- 
tance from the tip of the notch. The results are compared with stress values obtained 
by the boundary collocation technique reported by Gross (ref. 2) and are found to be in 
excellent agreement. Table I contains selected results of these stress computations. 

As expected, the stresses approach infinity near the tip of the notch. The square 
root singularity associated with the crack changes from 0. 500 to approximately 0. 488 
for a 60° notch angle as shown in table II. The order of stress singularity computed 
herein is compared with results given by Gross and Mendelson (ref. 8) and excellent 
agreement is obtained. 

The stress intensity factor Kj under mode I notch displacement is defined (ref. 2), 
in terms of the coordinate system shown in figure 2, as follows: 


K t = lim 
1 r — 0 



(19) 


From the known stress field in the vicinity of the tip of the notch the order of stress 
singularity n in equation (19) can be determined by plotting In o against In r, and 
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obtaining a least-squares fit of a straight line through the plotted points. The slope of 
the line gives n. 

Once the order of stress singularity is found, a plot of the function of r in equa- 
tion (19) as r — 0 yields the stress intensity factor K T . Figure 10 shows the variation 

A © 

of the dimensionless stress intensity factor Kj with notch depth for a 10 notch as 
compared to the analytical solution obtained by Gross and Mendelson (ref. 8). These 
analytical results generally give values higher by approximately 1 percent. 

Figure 11 shows the dimensionless elastic plane-stress displacements at the edge 
of the notch as a function of notch depth. The displacements were obtained by numerical 
integration of relation (22) along straight line paths. The comparison is made with 
available results by Gross (ref. 2) and are found to be in very good agreement. 

Finally, the relations between path independent Rice's integral J and the stress 
intensity factor Kj were checked. These relations were given by Rice (ref. 9) for a 
linear elastic body in the following form: 

2 

J = — — K? (for plane strain) (20a) 

E 1 


J = — Kj (for plane stress) (20b) 

E 1 

where J is defined as 



( 21 ) 


Here r is an arbitrary curve surrounding the notch tip. The integral is evaluated in a 
counterclockwise sense, starting from the notch surface. Strain energy density is de- 
fined as 


W(e ) = 
' ran' 


r mn 


and the traction vector is 
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while u A is a displacement vector defined in terms of strain tensor as 


e ii = “ ( u i i + u i i) 

2 1 »J L 1 


( 22 ) 


Equations (20) allow one to evaluate the stress intensity factors without detailed 
knowledge of the stress field very near the notch tip. 

For the plane- stress case, the comparison between J values obtained by numerical 
integration of relation (21) and by use of equation (20b) is shown in table in. 


CONCLUSIONS 

The results obtained for the elastic stress analysis of V-notched beams subjected 
to pure bending compared very well with the collocation method of solution. The bound- 
ary integral method appears to be well suited for the solution of problems with geomet- 
ric singularities. The dependence of the accuracy of the stress field in the vicinity of 
notch tip on the choice of nodal spacing has been demonstrated. The strategy of making 
this choice has been shown to be successful by comparison to results obtained by the 
collocation method of solution. The experience gained from those studies can now be 
utilized in performing the more complicated elastoplastic calculations. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 20, 1973, 

501-21. 
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APPENDIX - COEFFICIENTS OF THE STRESS EQUATIONS 


The coefficients appearing in stress equations (17) are given by the following rela- 
tions : 



ij 



In [(x. - |) 2 + (y i - t?) 2 ] 


2(Yi - nY 


(*i - i ) 2 + (Yi - V ) 2 


+ 1> ds 




2(x. - er 


(x. - ?r + (yj - tiY 


ds 




3x ( 


- (d..) = - 
2 1] 




2(x i - g) 2 

(Xi - I ) 2 + (Yi - r ,) 2 




+ lV ds 
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G-. =— - — (e..) = -8 I -2-J 
13 0x 9y ' 13 ' / 9n < 


(34 - iKyj - v) 


ds 


'i 


J(Xj - I) 2 + (yj - I )) 2 


J 3x 3y 


J h - «> 2 + - A 


ds 


V-E-,0.1 / 

1J Sx 9y 1J I 3n /„ >\2 . /__ _\2 


(Xj - ?) + (Yi - i?) 


ds 


K ii =- L ~ ( d ii) = - 2 
13 3x 3y 13 


(34 - 5)(y t - n) 

(34 - S) 2 + (y t - *?) 2 


ds 


The evaluation of these integrals is given in reference 5. 
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TABLE I. - DIMENSIONLESS ELASTIC x- DIRECTIONAL STRESSES a /a n AND y- DIRECTIONAL 

x 0 J 

STRESSES Vy/<>o ALONG x-AXIS (y = 0) IN THE VICINITY OF THE NOTCH FOR A 
SPECIMEN WITH A SINGLE EDGE NOTCH SUBJECTED TO PURE BENDING 
[Dimensionless load q = 1.0.] 

(a) Dimensionless notch depth a = 0.2 (b) Dimensionless notch depth a =• 0. 3 


x/a 

Notch angle, a , deg 



3 

10 

30 

60 

? 

X 

0.01 

7. 72 

7.88 

7.66 

6.72 

.02 

5. 45 

5.54 

5.33 

4. 80 

.04 

3. 80 

3. 85 

3.72 

3.40 

.06 

3. 06 

3.09 

3.00 

2.77 

. 10 

— 

2.30 

2.26 

2. 11 

.20 

— 

1.51 

1.50 

1.42 

? y 

0. 01 

7. 81 

7.99 

7.84 

7.81 

.02 

5.54 

5.65 

5.52 

5.54 

.04 

3. 88 

3.95 

3.86 

3. 89 

.06 

3. 13 

3. 18 

3. 10 

3. 14 

. 10 

— 

2.39 

2.33 

2.37 

.20 

— 

1.54 

1.51 

1.54 


x/a 

Notch angle, a, 

deg 


10 

30 

60 


r 

X 


0. 01 

7. 13 

7. 04 

6.28 

.02 

5.01 

4.93 

4. 47 

. 04 

3.47 

3. 43 

3. 15 

.06 

2. 77 

2. 75 

2. 56 

. 10 

2.07 

2. 06 

1.94 

.20 

1.35 

1.35 

1.29 


r, 

y 


0.01 

7.40 

7. 44 

7.36 

.02 

5.27 

5.27 

5.24 

. 04 

3. 73 

3. 72 

3.71 

.06 

3. 04 

3.03 

3. 03 

. 10 

2.33 

2. 27 

2.33 

.20 

1.61 

1. 60 

1.61 
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TABLE I. - Continued. DIMENSIONLESS ELASTIC x- DIRECTIONAL STRESSES a^/a Q AND y- DIRECTIONAL 
STRESSES CT y /a Q ALONG x-AXIS (y = 0) IN THE VICINITY OF THE NOTCH FOR A SPECIMEN WITH 
A SINGLE EDGE NOTCH SUBJECTED TO PURE BENDING 
[Dimensionless load q = 1.0. ] 

(c) Dimensionless notch depth a = 0. 4 (d) Dimensionless notch depth a = 0. 5 


x/a 

Notch angle a, deg 

10 

30 

60 


a 





K 


0.01 

8.82 

' 8.55 

7.53 

.02 

6.24 

6.07 

5.40 

. 04 

4.36 

4.26 

3.84 

.06 

3.51 

3.44 

3.13 

. 10 

2.63 

2.59 

j 2.39 

.20 

1.70 

1.68 

1.58 


a. 





y 


0.01 

8.74 

8.73 

8.66 

.02 

6. 13 

6.13 

6. 11 

.04 

4.23 

4.23 

4.25 

.06 

3.36 

3.36 

3.39 

. 10 

2.45 

2.45 

2.48 

.20 

1.44 

1.44 

i 

1.47 


x/a 

Notch angle, a, deg 

10 

30 

60 


a 





X 


0. 01 

10. 63 

10. 46 

8.98 

.02 

7.54 

7.42 

6.44 

.04 

5.28 

5.20 

4.58 

.06 

4.24 

4.19 

3.72 

.10 

3.15 

3. 13 

2.80 

.20 

1.96 

- 

1.94 

1.79 


a 





y 


0.01 

10.29 

10. 49 

10.21 

.02 

7.17 

7.29 

7.15 

.04 

4.86 

4.94 

4.88 

.06 

3.79 

3.85 

3.82 

.10 

2.63 

2.67 

2.68 

.20 

1.30 

1.32 

1.35 
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TABLE I. - Concluded. DIMENSIONLESS ELASTIC x- DIRECTIONAL 


STRESSES a x /cr 0 AND y- DIRECTIONAL STRESSES o / <Jq ALONG 
x-AXIS (y = 0) IN THE VICINITY OF THE NOTCH FOR A 
SPECIMEN WITH A SINGLE EDGE NOTCH 


SUBJECTED TO PURE BENDING 
(Dimensionless load q = 1.0J 
(e) Dimensionless notch depth a = 0. 6 


x/a 

Notch angle, a, deg 


10 

. . 

30 

60 


3 




X 


0.01 

13.80 

13.46 

11.59 

.02 

9. 78 

9.56 

8.30 

.04 

6. 80 

6.68 

5. 87 

.06 

5. 42 

5.33 

4.73 

. 10 

3.94 

3. 88 

3.50 

.20 

2. 27 

2.24 

2.04 


a 





y 


0.01 

13.06 

13.26 

12.98 

.02 

8.96 

9. 11 

8.99 

.04 

5.91 

6.01 

5.98 

.06 

4. 46 

4.53 

4.54 

. 10 

2.85 

2.91 

2.93 

.20 

.83 

.88 

.90 




TABLE II. - ORDER OF STRESS SINGULARITY n AT THE TIP OF THE NOTCH FOR SPECIMENS 


WITH A SINGLE EDGE NOTCH SUBJECTED TO PURE BENDING AND BEHAVING ELASTICALLY 


[Dimensionless load q = 1. o!] 


Dimension- 

Notch angle, a, deg 

less notch 













depth, 


3 



10 



30 



60 


a = — 






Source 

































Present 

Reference 8 

Present 

Reference 8 

Present 

Reference 8 

Present 

Reference 8 


report 



report 



report 



report 





Order of stress singularity at tip of notch, n 

0.2 


a 0. 5000 

0. 4990 

0. 4999 

0. 4986 

0. 4985 

0.4875 

0.4878 

.3 

0.4999 



.4999 



.4978 



.4896 



.4 




.5007 



.4989 



.4929 



.5 




.4999 



.4989 



.4886 



.6 




.4999 



.5010 



.4863 



.7 




.5010 



.5000 







a Value obtained for a = 0°. 


TABLE in. - DIMENSIONLESS ELASTIC PLANE-STRESS RICE'S 
INTEGRAL — FOR A SPECIMEN WITH A 10° EDGE 

fa 

NOTCH SUBJECTED TO PURE BENDING 


[Dimensionless load q = 1.0; Poisson's ratio fi = 0. 33 


Dimension- 

Dimensionless elastic plane 

-stress Rice's integral obtained 

less notch 
depth, 

By integration of 

From stress intensity 

~ a 

equation (21) 

factor Kj in equa- 

a = — 
w 


tion (20b) 

0.3 

1. 193 

1.212 

.5 

3.546 

3.408 


20 













= 1 ‘ 

subdivisions for P(x,y)cC. 


A 25 24 23 22 21 B 











OnW 



Figure 10. - Dimensionless stress intensity factor Kj for a specimen with a 10° 
notch subjected to pure bending and behaving elastically. Dimensionless load 
q = 1.0. 



Figure 11. - Dimensionless elastic plane-stress y-directional notch opening displacement 
for a specimen with a \(P edge notch subjected to pure bending. Dimensionless load 
q * 1. 0; Poisson's ratio p = 0. 33. 
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